Stability criterion for time-delay of cyber-physical power systems under distributed control

ABSTRACT

The present disclosure provides a stability criterion for time-delay of cyber-physical power systems under distributed control, which relates to a field of cyber-physical power systems technologies. The method first establishes an cyber side model of the cyber-physical power systems under distributed control and a physical power grid model of the cyber-physical power systems under distributed control respectively; then establishes simultaneous equations of the cyber side model and the physical power grid model to establish an unified differential algebraic equation model of the cyber-physical power systems under distributed control, so as to obtain a time-delay characteristic equation expression of the cyber-physical power systems under distributed control; determines a time-delay of each node, and solving the time-delay characteristic equation expression to obtain a maximum characteristic root, and performing the stability criterion for the time-delay of the cyber-physical power systems under distributed control according to a real part of the maximum characteristic root.

FIELD

The present disclosure relates to a field of cyber-physical power systems, and more particularly to a stability criterion for time-delay of cyber-physical power systems under distributed control.

BACKGROUND

Modern power systems have an in-depth interdependence between a cyber side and a physical side, and have become typical cyber-physical systems (CPS). The cyber-physical systems are complex systems constructing physical spaces and cyber spaces which are mapped each other and efficiently cooperated through advanced technologies such as perception, calculation, communication and control. In the field of smart grid technologies, the modern smart grid has been widely regarded as a kind of cyber-physical power systems due to physical power systems relying on many advanced information technology manners.

An influence of time-delays on the cyber-physical power systems cannot be ignored more and more. Dispatchers of the modern power systems have put forward extremely high requirements of time-delay on various optimized scheduling and control services. In particular, stability of the systems may be affected due to commands not being issued timely caused by the time-delay of the control services being too long. Therefore, there have been many methods to study a stability criterion of time-delay power systems under centralized control. With developments of distributed power sources, application prospects of various kinds of distributed architectures, such as distributed frequency control, voltage control in the cyber-physical power systems are becoming wider and wider. Therefore, the cyber-physical power systems under distributed control may be regarded as a kind of distributed cyber-physical power systems. Due to high dependence on the cyber systems, the stability of the distributed cyber-physical power systems has also become a key issue for distributed applications. Existing technical methods for a stability criterion for time-delay are mainly oriented to centralized cyber-physical power systems, whose time delays are relatively single (for example, Rekasius substitution is used to perform the stability criterion for the cyber-physical power systems under centralized/local time-delay power system stabilizers), but a criterion method for the distributed cyber-physical power systems is still blank. Therefore, the existing methods are not enough to analyze the stability for time-delay of the cyber-physical power systems under distributed control.

SUMMARY

The present disclosure seeks to solve at least one of the problems existing in the related art, and provides a stability criterion for time-delay of cyber-physical power systems under distributed control. The present disclosure is oriented to the stability for time-delay of the cyber-physical power systems under distributed control, and determines whether the cyber-physical power systems will be unstable under a specific time-delay. The method of the present disclosure may be applied to architecture design, communication networking and formulation of communication requirements of the power system under distributed control, and may clarify an influence of the distributed control on the stability of the cyber-physical power systems under different time-delays.

The present disclosure provides the stability criterion for time-delay of the cyber-physical power systems under distributed control. The method includes: establishing a cyber side model of the cyber-physical power systems under distributed control and a physical power grid model of the cyber-physical power systems under distributed control respectively; then establishing simultaneous equations of the cyber side model and the physical power grid model to establish an unified differential algebraic equation model of the cyber-physical power systems under distributed control, so as to obtain a time-delay characteristic equation expression of the cyber-physical power systems under distributed control; finally determining a time-delay of each node, and solving the time-delay characteristic equation expression to obtain a maximum characteristic root, and performing the stability criterion for time-delay of the cyber-physical power systems under distributed control according to a real part of the maximum characteristic root.

The method includes the following blocks.

1) The cyber side model of the cyber-physical power systems under distributed control is established. The detailed blocks are as follows.

1.1) An equation for a measuring link of each node of the cyber-physical power systems under distributed control is established.

The equation for the measuring link of a node i of the cyber-physical power systems under distributed control is as follows.

$\begin{matrix} {{z_{i}^{c}(t)} = {\begin{bmatrix} M_{x,i} & 0 \\ 0 & M_{y,i} \end{bmatrix}\begin{bmatrix} {x_{i}^{p}\left( {t - \tau_{m,i}} \right)} \\ {y_{i}^{p}\left( {t - \tau_{m,i}} \right)} \end{bmatrix}}} & (13) \end{matrix}$

where a superscript c denotes a cyber side of the cyber-physical power systems, a superscript p denotes a physical side of the cyber-physical power systems, x_(i) ^(p))(t−τ_(m,i)) denotes a physical state variable of the node i at time point t−τ_(m,i), y_(i) ^(p)(t−τ_(m,i)) denotes a physical algebraic variable of the node i at time point t−τ_(m,i), z_(i) ^(c)(t) denotes a cyber measuring variable of the node i at time point t, τ_(m,i) denotes a time-delay for measuring of the node i, M_(x,i) denotes an incidence matrix from the physical state variable of the node i to the cyber measuring variable of the node i, and M_(y,i) denotes an incidence matrix from the physical algebraic variable of the node i to the cyber measuring variable of the node i.

1.2) An equation for a communication iterative calculating link of each node of the cyber-physical power systems under distributed control is established.

The equation for the communication iterative calculating link of the node i of the cyber-physical power systems under distributed control is as follows.

{dot over (x)} _(i) ^(c)(t)=h _(i)(z _(i) ^(c)(t),x _(j∈N) _(i) ^(c)(t−τ _(c,ij)))  (14)

where x_(i) ^(c)(t) denotes a cyber state variable of the node i at time point t, {dot over (x)}_(i) ^(c)(t) denotes a differential of x_(i) ^(c)(t) with respect to time, N_(i) denotes a set of reference numbers of neighbor nodes communicating with the node i each other, x_(j∈N) _(i) ^(c)(t) denotes a cyber state variable of the neighbor node of the node i at time point t, τ_(c,ij), denotes a time-delay for communication between the node i and a neighbor node j of the node i, and h_(i) denotes an iterative correction function of the node i under distributed control.

1.3) An equation for a control link of each node of the cyber-physical power systems under distributed control is established.

The equation for the control link of the node i of the cyber-physical power systems under distributed control is as follows.

u _(i) ^(p)(t)=C _(i) ·x _(i) ^(c)(t−τ _(e,i))  (15)

where C_(i) denotes an incidence matrix from the cyber state variable of the node i to a physical control variable of the node i, u_(i) ^(p)(t) denotes the physical control variable of the node i, and τ_(e,i) denotes a time-delay for controlling of the node i.

1.4) Simultaneous equations of the equations for the measuring link, the equations for the communication iteration calculating link and the equations for the control link of all nodes together are established, and the cyber side model of the cyber-physical power systems under distributed control is established as follows.

$\begin{matrix} \left\{ \begin{matrix} {{\overset{.}{x}}_{i}^{c} = {h_{i}\left( {{z_{i}^{c}(t)},{x_{j \in N_{i}}^{c}\left( {t - \tau_{c,{ij}}} \right)}} \right)}} & {{i = 1},{2\mspace{14mu}\ldots}\mspace{14mu},N} \\ {0 = {u_{i}^{p} - {C_{i} \cdot {x_{i\;}^{c}\left( {t - \tau_{e,i}} \right)}}}} & {{i = 1},{2\mspace{14mu}\ldots}\mspace{14mu},N} \\ \begin{matrix} {0 = {{z_{i}^{c}(t)} - {{{diag}\left( {M_{x,i},M_{y,i}} \right)} \cdot}}} \\ \left\lbrack {{x_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T},{y_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T}} \right\rbrack^{T} \end{matrix} & \; \\ \; & {{i = 1},{2\mspace{14mu}\ldots}\mspace{14mu},N} \end{matrix} \right. & (16) \end{matrix}$

where diag( ) denotes a partitioned diagonal matrix, and N denotes the number of nodes for the distributed control.

2) The physical power grid model of the cyber-physical power systems under distributed control is established. The detailed blocks are as follows.

2.1) A physical side differential equation of the cyber-physical power systems under distributed control is established as follows.

{dot over (x)} ^(p)(t)=f ^(p)(x ^(p)(t),y ^(p)(t),u ^(p)(t))  (17)

where f^(p) denotes a dynamic equation of the power systems, x^(p)(t) denotes a vector of the physical state variables of the cyber-physical power systems at time point t and satisfies x^(p)(t)=[x₁ ^(pT)(t), x₂ ^(pT)(t), . . . , x_(N) ^(pT) (t)]^(T), a superscript T denotes transpose, {dot over (x)}^(p)(t) denotes a differential of x^(p)(t) with respect to time, y^(p)(t) denotes a vector of the physical algebraic variables at time point t, and u^(p)(t) denotes a vector of the physical control variables at time point t and satisfies u^(p)(t)=[u₁ ^(pT)(t), u₂ ^(pT)(t), . . . , u_(N) ^(pT)(t)]^(T).

2.2) A physical side algebraic equation of the cyber-physical power systems under distributed control is established as follows.

0=g ^(p)(x ^(p)(t),y ^(p)(t))  (18)

where g^(p) denotes a stability power flow equation of the power systems.

2.3) A physical side algebraic equation with time-delay of the cyber-physical power systems under distributed control is established as follows.

$\begin{matrix} \left\{ \begin{matrix} {0 = {g^{p}\left( {{x^{p}\left( {t - \tau_{m,i}} \right)},{y^{p}\left( {t - \tau_{m,i}} \right)}} \right)}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {0 = {g^{p}\left( {{x^{p}\left( {t - \tau_{e,i}} \right)},{y^{p}\left( {t - \tau_{e,i}} \right)}} \right)}} & {{i = 1},2,\ldots\mspace{14mu},N} \end{matrix} \right. & (19) \end{matrix}$

2.4) Simultaneous equations of formulas (5)-(7) to establishing the physical power grid model of the cyber-physical power systems under distributed control are established.

3) Simultaneous equations of the formula (4) of the cyber side model in the block 1) and the formulas (5)-(7) of the physical power grid model in the block 2) are established to obtain the unified differential algebraic equation model of the cyber-physical power systems under distributed control as follows.

$\begin{matrix} \left\{ \begin{matrix} {{{\overset{.}{x}}^{p}(t)} = {f^{p}\left( {{x^{p}(t)},{y^{p}(t)},{u^{p}(t)}} \right)}} & \; \\ {0 = {g^{p}\left( {{x^{p}(t)},{y^{p}(t)}} \right)}} & \; \\ {0 = {g^{p}\left( {{x^{p}\left( {t - \tau_{m,i}} \right)},{y^{p}\left( {t - \tau_{m,i}} \right)}} \right)}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {0 = {g^{p}\left( {{x^{p}\left( {t - \tau_{e,i}} \right)},{y^{p}\left( {t - \tau_{e,i}} \right)}} \right)}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {{{{\overset{.}{x}}_{i}^{c}(t)} = {h_{i}\left( {{z_{i}^{c}(t)},{\left( x_{\;^{j \in N_{i}}}^{c} \right)\left( {t - \tau_{c,{ij}}} \right)}} \right)}}\ } & {{i = 1},2,\ldots\mspace{14mu},N} \\ {0 = {{u_{i}^{p}(t)} - {C_{i} \cdot {x_{i}^{c}\left( {t - \tau_{e,i}} \right)}}}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ \begin{matrix} {0 = {{z_{i}^{c}(t)} - {{{diag}\left( {M_{x,i},M_{y,i}} \right)} \cdot}}} \\ \left\lbrack {{x_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T},{y_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T}} \right\rbrack^{T} \end{matrix} & \; \\ \; & {{i = 1},2,\ldots\mspace{14mu},N} \end{matrix} \right. & (20) \end{matrix}$

4) The time-delay characteristic equation expression of the cyber-physical power systems under distributed control is established. The detailed blocks are as follows.

4.1) An equilibrium point of the cyber-physical power systems under distributed control is calculated, and the formula (8) is solved by making the differentials {dot over (x)}^(p)(t) and {dot over (x)}_(i) ^(c)(t) in the formula (8) both being 0, to obtain a vector of the equilibrium points of the cyber-physical power systems at time point t, namely, x(t)=[x ^(pT)(t), x ₁ ^(cT)(t), . . . , x _(i) ^(cT)(t), . . . , x _(N) ^(cT)(t)]^(T), where x ^(pT)(t) denotes a vector of the equilibrium points of the physical state variables at time point t, x ^(cT) (t) denotes the equilibrium point of the cyber state variable of the node i at time point t.

4.2) Linearization expansion is performed on the formula (8) of the unified differential algebraic equation model of the cyber-physical power systems under distributed control at the equilibrium point, and the following linearized differential algebraic equations is obtained.

$\begin{matrix} \left\{ \begin{matrix} {{\Delta\;{{\overset{.}{x}}^{p}(t)}} = {{A_{0}\;\Delta\;{x^{p}(t)}} + {B_{0}\Delta\;{y^{p}(t)}} + {\sum\limits_{i = 1}^{N}{E_{i}\Delta\;{u_{i}^{p}(t)}}}}} & \; \\ {0 = {{C_{0}\Delta\;{x^{p}(t)}} + {D_{0}\Delta\;{y^{p}(t)}}}} & \; \\ {0 = {{C_{e,i}\Delta\;{x^{p}\left( {t - \tau_{e,i}} \right)}} + {D_{e,i}\Delta\;{y^{p}\left( {t - \tau_{e,i}} \right)}}}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {0 = {{C_{m,i}\Delta\;{x^{p}\left( {t - \tau_{m,i}} \right)}} + {D_{m,i}\Delta\;{y^{p}\left( {t - \tau_{m,i}} \right)}}}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {{\Delta\;{\overset{.}{x}}_{i}^{c}} = {{\sum\limits_{J \in N_{i}}{A_{ij}^{c}\Delta{x_{j}^{c}\left( {t - \tau_{c,{ij}}} \right)}}} + {B_{i}^{c}\Delta\;{z_{i}^{c}(t)}}}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {{\Delta\; u_{i}^{p}} = {C_{i}\Delta\;{x_{i}^{c}\left( {t - \tau_{e,i}} \right)}}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {{\Delta\; z_{i}^{c}} = {{K_{x,i}\  \cdot {x_{i}^{p}\left( {t - \tau_{m,i}} \right)}} + {K_{y,i} \cdot {y_{i}^{p}\left( {t - \tau_{m,i}} \right)}}}} & {{i = 1},2,\ldots\mspace{14mu},N} \end{matrix} \right. & (21) \end{matrix}$

where, Δ denotes a deviation value of the variable after performing the linearization expansion at the equilibrium point, A₀ denotes a Jacobi matrix obtained by taking derivative of f^(p) with respect to x^(p)(t), namely satisfies A₀=∂f^(p)/∂x^(p)(t), B₀ denotes a Jacobi matrix obtained by taking the derivative of f^(p) with respect to y^(p)(t), namely satisfies B₀=∂f^(p)/∂y^(p)(t), E_(i) denotes a Jacobi matrix obtained by taking the derivative of f^(p) with respect to u_(i) ^(p)(t), namely satisfies E_(i)=∂f^(p)/∂u_(i) ^(p)(t), i=1, . . . , N, C₀ denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to x^(p)(t), namely satisfies C₀=∂g^(p)/∂x^(p)(t), D₀ denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to y_(p)(t), namely satisfies D₀=∂g^(p)/∂y^(p)(t), C_(e,i) denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to x^(p)(t−τ_(e,i)), namely satisfies C_(e,i)=∂g^(p)/∂x^(p)(t−τ_(e,i)), D_(e,i) denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to y^(p)(t−τ_(e,i)), namely satisfies D_(e,i)=∂g^(p)/∂y^(p)(t−τ_(e,i)), C_(m,i) denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to x^(p) (t−τ_(m,i)), namely satisfies C_(m,i)=∂g^(p)/∂x^(p)(t−r_(m,i)), D_(m,i), denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to y^(p)(t−τ_(m,i)), namely satisfies D_(m,i)=∂g^(p)/∂y^(p)(t−τ_(m,i)), A_(ij) ^(c) denotes a Jacobi matrix obtained by taking the derivative of h_(i) with respect to x_(j) ^(c)(t−τ_(c,ij)), namely satisfies A_(ij) ^(c)=∂h_(i)/∂x_(j) ^(c)(t−τ_(c,ij)), B_(i) ^(c) denotes a Jacobi matrix obtained by taking the derivative of h_(i) with respect to z_(i) ^(c)(t), namely satisfies B_(i) ^(c)=∂h_(i)/∂z_(i) ^(c)(t), C_(i) remains unchanged, K_(x,i) denotes a Jacobi matrix obtained by taking the derivative of diag(M_(x,i),M_(y,i)) with respect to x_(i) ^(p)(t−τ_(m,i)), namely satisfies K_(x,i)=∂diag(M_(x,i),M_(y,i))/∂x_(i) ^(p)(t−τ_(m,i)), and K_(y,i) denotes a Jacobi matrix obtained by taking the derivative of diag(M_(x,i),M_(y,i)) with respect to y_(i) ^(p)(t−τ_(m,i)), namely satisfies K_(y,i)=∂diag(M_(x,i),M_(y,i))/∂y_(i) ^(p)(t−τ_(m,i)).

4.3) Algebraic equations in the formula (9) of the linearized differential algebraic equations are eliminated, and the formula (9) is equivalently transformed into time-delay differential equations of the cyber-physical power systems under distributed control as follows:

$\begin{matrix} \left\{ \begin{matrix} {{\Delta\;\overset{.}{x}\;(t)} = {{\begin{bmatrix} {A_{0} - {B_{0}D_{0}^{- 1}C_{0}}} & 0 \\ 0 & 0 \end{bmatrix}\ \Delta\;{x(t)}} + {\sum\limits_{i = 1}^{N}\ {{\begin{bmatrix} 0 & {E_{i}C_{i}} \\ 0 & 0 \end{bmatrix} \cdot \Delta}\;{x\left( {t - \tau_{e,i}} \right)}}} +}} \\ {{\sum\limits_{i = 1}^{N}{\sum\limits_{j \in N_{i}}{{\begin{bmatrix} 0 & 0 \\ 0 & A_{ij}^{c} \end{bmatrix} \cdot \Delta}\;{x\left( {t - \tau_{c,{ij}}} \right)}}}} +} \\ {\sum\limits_{i = 1}^{N}\ {{\begin{bmatrix} 0 & 0 \\ {{B_{i}^{c}K_{x,i}} - {B_{i}^{c}K_{y,i}D_{m,i}^{- 1}C_{m,i}}} & 0 \end{bmatrix} \cdot \Delta}\;{x\left( {t - \tau_{m,i}} \right)}}} \end{matrix} \right. & (22) \end{matrix}$

where a superscript ⁻¹ of a matrix denotes an inverse matrix of this matrix, x denotes a composite vector of the physical state variables and the cyber state variables, and satisfies x=[x^(pT), x₁ ^(cT), . . . , x_(i) ^(cT), . . . x_(N) ^(cT)]^(T).

4.4) The following time-delay characteristic equation of the cyber-physical power systems under distributed control is established.

$\begin{matrix} {{\det\begin{pmatrix} {{s \cdot I} - {\overset{\sim}{A}}_{0} - {\sum\limits_{i = 1}^{N_{c}}{T_{m,i} \cdot e^{{- s}\;\tau_{m,i}}}} -} \\ {{\sum\limits_{i = 1}^{N_{c}}{\sum\limits_{j \in N_{i}}{T_{ij} \cdot e^{{- s}\;\tau_{c,{ij}}}}}} - {\sum\limits_{i = 1}^{N_{c}}{T_{e,i} \cdot e^{{- s}\;\tau_{e,i}}}}} \end{pmatrix}} = 0} & (23) \end{matrix}$

where det denotes a determinant, s denotes a characteristic root, e denotes a natural exponential, and Ã₀, T_(ij), T_(m,i), T_(e,i) respectively satisfy:

$\begin{matrix} \left\{ \begin{matrix} {{{\overset{\sim}{A}}_{0} = \begin{bmatrix} {A_{0} - {B_{0}D_{0}^{- 1}C_{0}}} & 0 \\ 0 & 0 \end{bmatrix}},{T_{ij} = \begin{bmatrix} 0 & 0 \\ 0 & A_{ij}^{c} \end{bmatrix}}} \\ {{T_{m,i} = \begin{bmatrix} 0 & 0 \\ {{B_{i}^{c}K_{x,i}} - {B_{i}^{c}K_{y,i}D_{m,i}^{- 1}C_{m,i}}} & 0 \end{bmatrix}},{T_{e,i} = \begin{bmatrix} 0 & {E_{i}c_{i}} \\ 0 & 0 \end{bmatrix}}} \end{matrix} \right. & (24) \end{matrix}$

5) The time-delay of each node is determined and the stability criterion for the time-delay of the cyber-physical power systems under distributed control is performed. The detailed blocks are as follows.

5.1) The time-delay for measuring of each node τ_(m,i), i=1, 2 . . . , N, the time-delay for controlling of each node τ_(e,i)=1, 2 . . . , N and the time-delay τ_(c,ij) for communication between each node i=1, 2, . . . , N and neighbor nodes j∈N, of the node are set, then the formula (11) of the time-delay characteristic equation is solved to obtain all characteristic roots of the equation. A characteristic root having the maximum real part is denoted as λ₁, and the real part of λ₁ is Re(λ₁).

5.2) It is determined whether Re(λ₁) is greater than 0. In response to Re(λ₁) being greater than 0, the cyber-physical power systems under distributed control is determined to be unstable. In response to Re(λ₁) being less than or equal to 0, the cyber-physical power systems under distributed control is determined to be stable. The stability criterion for the time-delay completes.

Characteristics and beneficial effects of the present disclosure are:

The present disclosure provides a stability criterion for time-delay of cyber-physical power systems under distributed control. The method of the present disclosure may be applied to the architecture design, the communication networking and the formulation of the communication requirements of the power system under distributed control, and may clarify the influence of the distributed control on the stability of the cyber-physical power systems under different time-delays.

DETAILED DESCRIPTION

The present disclosure provides a stability criterion for time-delay of cyber-physical power systems under distributed control. The following will further explain in combination with specific embodiments.

The present disclosure provides a stability criterion for time-delay of cyber-physical power systems under distributed control. First a cyber side model of the cyber-physical power systems under distributed control is established. The cyber-side model includes an equation for a measuring link, an equation for a communication iterative calculating link and an equation for a control link of the cyber-physical power systems under distributed control. A physical power grid model of the cyber-physical power systems under distributed control is established. The model includes a differential equation reflecting a dynamic process of a power grid and an algebraic equation reflecting a connection relationship of the power grid and a steady state process. Then simultaneous equations of the cyber side model and the physical power grid model is established to establish a unified differential algebraic equation model of the cyber-physical power systems under distributed control, so as to obtain a time-delay characteristic equation expression of the cyber-physical power systems under distributed control. Finally, a time-delay of each node is determined, and the time-delay characteristic equation expression is solved to obtain a maximum characteristic root, and the stability criterion for the time-delay of the cyber-physical power systems under distributed control is performed according to a real part of the maximum characteristic root. The method of the present disclosure may be applied to architecture design, communication networking and formulation of communication requirements of the power system under distributed control, and may clarify an influence of the distributed control on the stability of the cyber-physical power systems under different time-delays.

The method includes the following blocks.

1) The cyber side model of the cyber-physical power systems under distributed control is established. The cyber side model includes the equation for the measuring link, the equation for the communication iterative calculating link, and the equation for the control link of the cyber-physical power systems under distributed control. The detailed blocks are as follows.

1.1) An equation for a measuring link of each node of the cyber-physical power systems under distributed control is established.

The equation for the measuring link of a node i of the cyber-physical power systems under distributed control is as follows.

$\begin{matrix} {{z_{i}^{c}(t)} = {\begin{bmatrix} M_{x,i} & 0 \\ 0 & M_{y,i} \end{bmatrix} \cdot \begin{bmatrix} {x_{i}^{p}\left( {t - \tau_{m,i}} \right)} \\ {y_{i}^{p}\left( {t - \tau_{m,i}} \right)} \end{bmatrix}}} & (25) \end{matrix}$

where a superscript c denotes a cyber side of the cyber-physical power systems, a superscript p denotes a physical side of the cyber-physical power systems, x_(i) ^(p)(t−τ_(m,i)) denotes a physical state variable of the node i at time point t−τ_(m,i) (for example, phase angle of the node, frequency, voltage, etc., the physical state variable may be selected one or more from these), y_(i) ^(p)(t−τ_(m,i)) denotes a physical algebraic variable of the node i at time point t−τ_(m,i) (For example, active power output of the node, reactive power output, active power load, reactive power load, etc., the physical algebraic variable may be selected one or more from these), z_(i) ^(c)(t) denotes a cyber measuring variable of the node i at time point t (the physical state variable and physical algebraic variable of the node i are converted into digital quantities on the cyber side through measuring), τ_(m,i) denotes a time-delay for measuring of the node i, M_(x,i) denotes an incidence matrix from the physical state variable of the node i to the cyber measuring variable of the node i, where x denotes the corresponding physical state variable, i denotes the reference number corresponding to the node, when a fault is not considered, M_(x,i) denotes a diagonal matrix with diagonal elements all 1 (the number of rows of the matrix is the same as a dimension of z_(i) ^(c)(t), and the number of columns of the matrix is the same as a dimension of x_(i) ^(p)(t−τ_(m,i))), M_(y,i) denotes an incidence matrix from the physical algebraic variable of the node i to the cyber measuring variable of the node i, where y denotes the corresponding physical algebraic variable, i denotes the reference number corresponding to the node, when the fault is not considered, M_(y,i) denotes a diagonal matrix with diagonal elements all 1 (the number of rows of the matrix is the same as a dimension of z_(i) ^(c)(t), and the number of columns of the matrix is the same as a dimension of y_(i) ^(p)(t−τ_(m,i))).

1.2) An equation for a communication iterative calculating link of each node of the cyber-physical power systems under distributed control is established.

The equation for the communication iterative calculating link of the node i of the cyber-physical power systems under distributed control is as follows.

{dot over (x)} _(i) ^(c)(t)=h _(i)(z _(i) ^(c)(t),x _(j∈N) _(i) ^(c)(t−τ _(c,ij)))  (26)

where x_(i) ^(c)(t) denotes a cyber state variable of the node i at time point t (for example, a tiny increasing rate or an output of a generator of the node i for advanced application decision-making on the cyber side depends on parameters of the specific algorithm iteration under distributed control), {dot over (x)}_(i) ^(c)(t) denotes a differential of x_(i) ^(c)(t) with respect to time, namely {dot over (x)}_(i) ^(c)(t)=dx_(i) ^(c)(t)/dt, N_(i) denotes a set of reference numbers of neighbor nodes communicating with the node i each other, x_(j∈N) _(i) ^(c)(t) denotes a cyber state variable of the neighbor node of the node i at time point t, τ_(c,ij) denotes a time-delay for communication between the node i and a neighbor node j of the node i, and h_(i) denotes an iterative correction function of the node i under distributed control.

1.3) An equation for a control link of each node of the cyber-physical power systems under distributed control is established.

The equation for the control link of the node i of the cyber-physical power systems under distributed control is as follows.

u _(i) ^(p)(t)=C _(i) ·x _(i) ^(c)(t−τ _(e,i))  (27)

where C_(i) denotes an incidence matrix from the cyber state variable of the node i to a physical control variable of the node i (the number of rows of the matrix is the same as a dimension of u_(i) ^(p)(t), and the number of columns of the matrix is the same as a dimension of x_(i) ^(c)(t−τ_(e,i))), when the fault is not considered, C_(i) denotes a diagonal matrix with diagonal elements all 1, u_(i) ^(p)(t) denotes the physical control variable of the node i (corresponding to an adjustable device in the physical power grid, such as the active power output of the generator, etc.), and τ_(e,i) denotes a time-delay for controlling of the node i.

1.4) Simultaneous equations of the formula (1) to (3) of the equations for the measuring link, the equations for the communication iteration calculating link and the equations for the control link of all nodes (i=1, . . . , N, N denotes the number of nodes for the distributed control) together are established, and the cyber side model of the cyber-physical power systems under distributed control is established as follows.

$\begin{matrix} \left\{ \begin{matrix} {{\overset{.}{x}}_{i}^{c} = {h_{i}\ \left( {{z_{i}^{c}(t)},{x_{\;^{j \in N_{i}}}^{c}\ \left( {t - \tau_{c,{ij}}} \right)}} \right)}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {0 = {u_{i}^{p} - {C_{i} \cdot {x_{i}^{c}\left( {t - \tau_{e,i}} \right)}}}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ \begin{matrix} {0 = {{z_{i}^{c}(t)} - {{{diag}\left( {M_{x,i},M_{y,i}} \right)} \cdot}}} \\ \left\lbrack {{x_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T},{y_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T}} \right\rbrack^{T} \end{matrix} & \; \\ \; & {{i = 1},2,\ldots\mspace{14mu},N} \end{matrix} \right. & (28) \end{matrix}$

where diag( ) denotes a partitioned diagonal matrix.

2) The physical power grid model of the cyber-physical power systems under distributed control is established. The model includes the differential equation reflecting the dynamic process of the power grid and the algebraic equation reflecting the connection relationship of the power grid and the steady state process. The detailed blocks are as follows.

2.1) A physical side differential equation of the cyber-physical power systems under distributed control is established as follows.

{dot over (x)} ^(p)(t)=f ^(p)(x ^(p)(t),y ^(p)(t),u ^(p)(t))  (29)

where f^(p) denotes a dynamic equation of the power systems, x^(p)(t) denotes a vector of the physical state variables of the cyber-physical power systems at time point t and satisfies x^(p)(t)=[x₁ ^(pT)(t), x₂ ^(pT)(t), . . . , x_(N) ^(pT) (t)]^(T), a superscript T denotes transpose, {dot over (x)}^(p)(t) denotes a differential of x^(p)(t) with respect to time, y^(p) (t) denotes a vector of the physical algebraic variables at time point t, u^(p)(t) denotes a vector of the physical control variables at time point t and satisfies u^(p)(t)=[u₁ ^(pT)(t), u₂ ^(pT)(t), . . . , u_(N) ^(pT)(t)]^(T).

2.2) A physical side algebraic equation of the cyber-physical power systems under distributed control is established as follows.

0=g ^(p)(x ^(p)(t),y ^(p)(t))  (30)

where g^(p) denotes a stability power flow equation of the power systems.

2.3) A physical side algebraic equation with time-delay of the cyber-physical power systems under distributed control is established as follows.

$\begin{matrix} \left\{ \begin{matrix} {0 = {g^{p}\left( {{x^{p}\left( {t - \tau_{m,i}} \right)},{y^{p}\left( {t - \tau_{m,i}} \right)}} \right)}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {0 = {g^{p}\left( {{x^{p}\left( {t - \tau_{e,i}} \right)},{y^{p}\left( {t - \tau_{e,i}} \right)}} \right)}} & {{i = 1},2,\ldots\mspace{14mu},N} \end{matrix} \right. & (31) \end{matrix}$

2.4) Simultaneous equations of formulas (5)-(7) to establishing the physical power grid model of the cyber-physical power systems under distributed control are established.

3) The unified differential algebraic equation model of the cyber-physical power systems under distributed control is established.

Simultaneous equations of the formula (4) of the cyber side model in the block 1) and the formulas (5)-(7) of the physical power grid model in the block 2) are established, to obtain the unified differential algebraic equation model of the cyber-physical power systems under distributed control as follows.

$\begin{matrix} \left\{ \begin{matrix} {{{\overset{.}{x}}^{p}(t)} = {f^{p}\left( {{x^{p}(t)},{y^{p}(t)},{u^{p}(t)}} \right)}} & \; \\ {0 = {g^{p}\left( {{x^{p}(t)},{y^{p}(t)}} \right)}} & \; \\ {0 = {g^{p}\left( {{x^{p}\left( {t - \tau_{m,i}} \right)},{y^{p}\left( {t - \tau_{m,i}} \right)}} \right)}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {0 = {g^{p}\left( {{x^{p}\left( {t - \tau_{e,i}} \right)},{y^{p}\left( {t - \tau_{e,i}} \right)}} \right)}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {{{\overset{.}{x}}_{i}^{c}(t)} = {h_{i}\left( {{z_{i}^{c}(t)},{x_{\;^{j \in N_{i}}}^{c}\ \left( {t - \tau_{c,{ij}}} \right)}} \right)}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {0 = {{u_{i}^{p}(t)} - {C_{i} \cdot {x_{i}^{c}\left( {t - \tau_{e,i}} \right)}}}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ \begin{matrix} {0 = {{z_{i}^{c}(t)} - {{{diag}\left( {M_{x,i},M_{y,i}} \right)} \cdot}}} \\ \left\lbrack {{x_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T},{y_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T}} \right\rbrack^{T} \end{matrix} & \; \\ \; & {{i = 1},2,\ldots\mspace{14mu},N} \end{matrix} \right. & (32) \end{matrix}$

4) The time-delay characteristic equation expression of the cyber-physical power systems under distributed control is established. The detailed blocks are as follows.

4.1) An equilibrium point of the cyber-physical power systems under distributed control is calculated, and the formula (8) is solved by making the differentials {dot over (x)}^(p)(t) and {dot over (x)}_(i) ^(c)(t) in the formula (8) both being 0, to obtain a vector of the equilibrium points (that is, a point at which the system may maintain stability without disturbance) of the cyber-physical power systems at time point t, namely, x(t)=[x ^(pT)(t), x ₁ ^(cT)(t), . . . , x _(i) ^(cT)(t), . . . , x _(N) ^(cT)(t)]^(T), where x ^(pT)(t) denotes a vector of the equilibrium points of the physical state variables at time point t (the vector includes the equilibrium point of each physical state variable), x _(i) ^(cT)(t) denotes the equilibrium point of the information state variable of the node i at time point t.

4.2) Linearization expansion is performed on the formula (8) of the unified differential algebraic equation model of the cyber-physical power systems under distributed control at the equilibrium point, and the following linearized differential algebraic equations may be obtained.

$\begin{matrix} \left\{ \begin{matrix} {{\Delta\;{{\overset{.}{x}}^{p}(t)}} = {{A_{0}\;\Delta\;{x^{p}(t)}} + {B_{0}\Delta\;{y^{p}(t)}} + {\sum\limits_{i = 1}^{N}{E_{i}\Delta\;{u_{i}^{p}(t)}}}}} & \; \\ {0 = {{C_{0}\Delta\;{x^{p}(t)}} + {D_{0}\Delta\;{y^{p}(t)}}}} & \; \\ {0 = {{C_{e,i}\Delta\;{x^{p}\left( {t - \tau_{e,i}} \right)}} + {D_{e,i}\Delta\;{y^{p}\left( {t - \tau_{e,i}} \right)}}}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {0 = {{C_{m,i}\Delta\;{x^{p}\left( {t - \tau_{m,i}} \right)}} + {D_{m,i}\Delta\;{y^{p}\left( {t - \tau_{m,i}} \right)}}}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {{\Delta\;{\overset{.}{x}}_{i}^{c}} = {{\sum\limits_{J \in N_{i}}{A_{ij}^{c}\Delta{x_{j}^{c}\left( {t - \tau_{c,{ij}}} \right)}}} + {B_{i}^{c}\Delta\;{z_{i}^{c}(t)}}}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {{\Delta\; u_{i}^{p}} = {C_{i}\Delta\;{x_{i}^{c}\left( {t - \tau_{e,i}} \right)}}} & {{i = 1},2,\ldots\mspace{14mu},N} \\ {{\Delta\; z_{i}^{c}} = {{K_{x,i}\  \cdot {x_{i}^{p}\left( {t - \tau_{m,i}} \right)}} + {K_{y,i} \cdot {y_{i}^{p}\left( {t - \tau_{m,i}} \right)}}}} & {{i = 1},2,\ldots\mspace{14mu},N} \end{matrix} \right. & (33) \end{matrix}$

where Δ denotes a deviation value of the variable after performing the linearization expansion at the equilibrium point, A₀ denotes a Jacobi matrix obtained by taking derivative of f^(p) with respect to x^(p)(t), namely satisfies A₀=∂f^(p)/∂x^(p)(t), B₀ denotes a Jacobi matrix obtained by taking the derivative of f^(p) with respect to y^(p)(t), namely satisfies B₀=∂f^(p)/∂y^(p)(t), E_(i) denotes a Jacobi matrix obtained by taking the derivative of f^(p) with respect to u_(i) ^(p)(t), namely satisfies E_(i)=∂f^(p)/∂u_(i) ^(p)(t), i=1, . . . , N, C₀ denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to x^(p)(t), namely satisfies C₀=∂g^(p)/∂x^(p)(t), D₀ denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to y_(p)(t), namely satisfies D₀=∂g^(p)/∂y^(p)(t), C_(e,i) denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to x^(p)(t−τ_(e,i)), namely satisfies C_(e,i)=∂g^(p)/∂x^(p)(t−τ_(e,i)), D_(e,i) denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to y^(p)(t−τ_(e,i)), namely satisfies D_(e,i)=∂g^(p)/∂y^(p)(t−τ_(e,i)), C_(m,i) denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to x^(p) (t−τ_(m,i)), namely satisfies C_(m,i)=∂g^(p)/∂x^(p)(t−r_(m,i)), D_(m,i), denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to y^(p)(t−τ_(m,i)), namely satisfies D_(m,i)=∂g^(p)/∂y^(p)(t−τ_(m,i)), A_(ij) ^(c) denotes a Jacobi matrix obtained by taking the derivative of h_(i) with respect to x_(j) ^(c)(t−τ_(c,ij)), namely satisfies A_(ij) ^(c)=∂h_(i)/∂x_(j) ^(c)(t−τ_(c,ij)), B_(i) ^(c) denotes a Jacobi matrix obtained by taking the derivative of h_(i) with respect to z_(i) ^(c)(t), namely satisfies B_(i) ^(c)=∂h_(i)/∂z_(i) ^(c)(t), C_(i) remains unchanged, K_(x,i) denotes a Jacobi matrix obtained by taking the derivative of diag(M_(x,i),M_(y,i)) with respect to x_(i) ^(p)(t−τ_(m,i)), namely satisfies K_(x,i)=∂diag(M_(x,i),M_(y,i))/∂x_(i) ^(p)(t−τ_(m,i)), and K_(y,i) denotes a Jacobi matrix obtained by taking the derivative of diag(M_(x,i),M_(y,i)) with respect to y_(i) ^(p)(t−τ_(m,i)), namely satisfies K_(y,i)=∂diag(M_(x,i),M_(y,i))/∂y_(i) ^(p)(t−τ_(m,i)).

4.3) Algebraic equations in the formula (9) of the linearized differential algebraic equations are eliminated, and the formula (9) is equivalently transformed into time-delay differential equations of the cyber-physical power systems under distributed control as follows.

$\begin{matrix} \left\{ \begin{matrix} {{\Delta\;\overset{.}{x}\;(t)} = {{\begin{bmatrix} {A_{0} - {B_{0}D_{0}^{- 1}C_{0}}} & 0 \\ 0 & 0 \end{bmatrix}\ \Delta\;{x(t)}} + {\sum\limits_{i = 1}^{N}\ {{\begin{bmatrix} 0 & {E_{i}C_{i}} \\ 0 & 0 \end{bmatrix} \cdot \Delta}\;{x\left( {t - \tau_{e,i}} \right)}}} +}} \\ {{\sum\limits_{i = 1}^{N}{\sum\limits_{j \in N_{i}}{{\begin{bmatrix} 0 & 0 \\ 0 & A_{ij}^{c} \end{bmatrix} \cdot \Delta}\;{x\left( {t - \tau_{c,{ij}}} \right)}}}} +} \\ {\sum\limits_{i = 1}^{N}\ {{\begin{bmatrix} 0 & 0 \\ {{B_{i}^{c}K_{x,i}} - {B_{i}^{c}K_{y,i}D_{m,i}^{- 1}C_{m,i}}} & 0 \end{bmatrix} \cdot \Delta}\;{x\left( {t - \tau_{m,i}} \right)}}} \end{matrix} \right. & (34) \end{matrix}$

where a superscript ⁻¹ of a matrix denotes an inverse matrix of this matrix, x denotes a composite vector of the physical state variables and the information state variables, and satisfies x=[x^(pT), x₁ ^(cT), . . . , x_(i) ^(cT), . . . x_(N) ^(cT)]^(T).

4.4) The following time-delay characteristic equation of the cyber-physical power systems under distributed control is established.

$\begin{matrix} {{\det\begin{pmatrix} {{s \cdot I} - {\overset{\sim}{A}}_{0} - {\sum\limits_{i = 1}^{N_{c}}{T_{m,i} \cdot e^{{- s}\;\tau_{m,i}}}} -} \\ {{\sum\limits_{i = 1}^{N_{c}}{\sum\limits_{j \in N_{i}}{T_{ij} \cdot e^{{- s}\;\tau_{c,{ij}}}}}} - {\sum\limits_{i = 1}^{N_{c}}{T_{e,i} \cdot e^{{- s}\;\tau_{e,i}}}}} \end{pmatrix}} = 0} & (35) \end{matrix}$

where det denotes a determinant, s denotes a characteristic root, e denotes a natural exponential, and Ã₀, T_(ij), T_(m,i), T_(e,i) respectively satisfy:

$\begin{matrix} \left\{ \begin{matrix} {{{\overset{\sim}{A}}_{0} = \begin{bmatrix} {A_{0} - {B_{0}D_{0}^{- 1}C_{0}}} & 0 \\ 0 & 0 \end{bmatrix}},{T_{ij} = \begin{bmatrix} 0 & 0 \\ 0 & A_{ij}^{c} \end{bmatrix}}} \\ {{T_{m,i} = \begin{bmatrix} 0 & 0 \\ {{B_{i}^{c}K_{x,i}} - {B_{i}^{c}K_{y,i}D_{m,i}^{- 1}C_{m,i}}} & 0 \end{bmatrix}},{T_{e,i} = \begin{bmatrix} 0 & {E_{i}c_{i}} \\ 0 & 0 \end{bmatrix}}} \end{matrix} \right. & (36) \end{matrix}$

5) The time-delay of each node is determined and the stability criterion for the time-delay of the cyber-physical power systems under distributed control is performed. The detailed blocks are as follows.

5.1) The time-delay for measuring of each node τ_(m,i), i=1, 2 . . . , N, the time-delay for controlling of each node τ_(e,i), i=1, 2 . . . , N and the time-delay τ_(c,ij) for communication between each node i=1, 2, . . . , N and neighbor nodes j∈N_(i) of the node are set, then the formula (11) of the time-delay characteristic equation is solved, to obtain all characteristic roots of the equation. A characteristic root having the maximum real part is denoted as λ₁, and the real part of λ₁ is Re(λ₁).

5.2) It is determined whether Re(λ₁) is greater than 0. In response to Re(λ₁) being greater than 0, the cyber-physical power systems under distributed control is determined to be unstable. In response to Re(λ₁) being less than or equal to 0, the cyber-physical power systems under distributed control is determined to be stable. The stability criterion for the time-delay completes.

The following further describes the present disclosure in detail in combination with an embodiment as follows.

This embodiment provides the stability criterion for time-delay of cyber-physical power systems under distributed control. The method includes the following blocks.

1) The cyber side model of the cyber-physical power systems under distributed control is established. The cyber side model includes the equation for the measuring link, the equation for the communication iterative calculating link, and the equation for the control link of the cyber-physical power systems under distributed control. The detailed blocks are as follows.

1.1) The equation for the measuring link of each node of the cyber-physical power systems under distributed control is established.

The equation for the measuring link of the node i of the cyber-physical power systems under distributed control is as follows.

$\begin{matrix} {{z_{i}^{c}(t)} = {\begin{bmatrix} M_{x,i} & 0 \\ 0 & M_{y,i} \end{bmatrix} \cdot \begin{bmatrix} {x_{i}^{p}\left( {t - \tau_{m,i}} \right)} \\ {y_{i}^{p}\left( {t - \tau_{m,i}} \right)} \end{bmatrix}}} & (37) \end{matrix}$

where the superscript c denotes an cyber side of the cyber-physical power systems, the superscript p denotes the physical side of the cyber-physical power systems, x_(i) ^(p)(t−τ_(m,i)) denotes the physical state variable of the node i at time point t−τ_(m,i) (in this embodiment the phase angle of the node and the frequency of the node in a second-order equation of a dynamic power system are selected), y_(i) ^(p)(t−τ_(m,i)) denotes the physical algebraic variable of the node i at time point t−τ_(m,i) (in this embodiment the active power output of the node, the active power load of the node are selected), z_(i) ^(c)(t) denotes an information measuring variable of the node i at time point t (the physical state variable and physical algebraic variable of the node i are converted into digital quantities on the cyber side through measuring), τ_(m,i) denotes the time-delay for measuring of the node i, M_(x,i) denotes an incidence matrix from the physical state variable of the node i to the information measuring variable of the node i, where x denotes the corresponding physical state variable, i denotes the reference number corresponding to the node, when the fault is not considered, M_(x,i) denotes the diagonal matrix with diagonal elements all 1 (the number of rows of the matrix is the same as the dimension of z_(i) ^(c)(t), and the number of columns of the matrix is the same as the dimension of x_(i) ^(p)(t−τ_(m,i))), M_(y,i) denotes an incidence matrix from the physical algebraic variable of the node i to the information measuring variable of the node i, where y denotes the corresponding physical algebraic variable, i denotes the reference number corresponding to the node, when the fault is not considered, M_(y,i) denotes the diagonal matrix with diagonal elements all 1 (the number of rows of the matrix is the same as the dimension of z_(i) ^(c)(t), and the number of columns of the matrix is the same as the dimension of y_(i) ^(p)(t−τ_(m,i))).

1.2) The equation for the communication iterative calculating link of each node of the cyber-physical power systems under distributed control is established.

The equation for the communication iterative calculating link of the node i of the cyber-physical power systems under distributed control is as follows.

{dot over (x)} _(i) ^(c)(t)=h _(i)(z _(i) ^(c)(t),x _(j∈N) _(i) ^(c)(t−τ _(c,ij)))  (38)

where x_(i) ^(c)(t) denotes an information state variable of the node i at time point t (in this embodiment the tiny increasing rate of the generator of is selected to realize economic distribution of the second frequency modulation stage according to distributed frequency control), {dot over (x)}_(i) ^(c)(t) denotes the differential of x_(i) ^(c)(t) with respect to time, namely {dot over (x)}_(i) ^(c)(t)=dx_(i) ^(c)(t)/dt, N_(i) denotes the set of reference numbers of neighbor nodes communicating with the node i each other, x_(j∈N) _(i) ^(c)(t) denotes an information state variable of the neighbor node of the node i at time point t, τ_(c,ij) denotes the time-delay for communication between the node i and the neighbor node j of the node i, and h_(i) denotes an iterative correction function of the node i under distributed control.

1.3) The equation for the control link of each node of the cyber-physical power systems under distributed control is established.

The equation for the control link of the node i of the cyber-physical power systems under distributed control is as follows.

u _(i) ^(p)(t)=C _(i) ·x _(i) ^(c)(t−τ _(e,i))  (39)

where C_(i) denotes an incidence matrix from the information state variable of the node i to the physical control variable of the node i (the number of rows of the matrix is the same as the dimension of u_(i) ^(p)(t), and the number of columns of the matrix is the same as the dimension of x_(i) ^(c)(t−τ_(e,i))), when the fault is not considered, C_(i) denotes the diagonal matrix with diagonal elements all 1, u_(i) ^(p)(t) denotes the physical control variable of the node i (corresponding to an adjustable device in the physical power grid, such as the active power output of the generator, etc.), and τ_(e,i) denotes the time-delay for controlling of the node i.

1.4) The simultaneous equations of the formula (1) to (3) of the equations for the measuring link, the equations for the communication iteration calculating link and the equations for the control link of all nodes (i=1, . . . , N, N denotes the number of nodes for the distributed control) together are established, and the cyber side model of the cyber-physical power systems under distributed control is established as follows.

$\begin{matrix} \left\{ \begin{matrix} {{{\overset{\cdot}{x_{i}}}^{c} = {{{h_{i}\ \left( {{z_{i}^{c}(t)},\ {x_{j \in N_{i}}^{c}\ \left( {t - \tau_{c,{ij}}} \right)}} \right)}\ i} = 1}},{2\mspace{14mu}\ldots}\mspace{11mu},N} \\ {{0 = {{u_{i}^{p} - {{C_{i} \cdot {x_{i}^{c}\left( {t - \tau_{e,i}} \right)}}\ i}} = 1}},{2\mspace{14mu}\ldots}\mspace{11mu},N} \\ {0 = {{z_{i}^{c}(t)} - {dia{{g\left( {M_{x,i},M_{y,i}} \right)} \cdot \left\lbrack {{x_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T},{y_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T}} \right\rbrack^{T}}}}} \\ {{i = 1},{2\mspace{14mu}\ldots}\mspace{11mu},N} \end{matrix} \right. & (40) \end{matrix}$

where diag( ) denotes the partitioned diagonal matrix.

2) The physical power grid model of the cyber-physical power systems under distributed control is established. The model includes the differential equation reflecting the dynamic process of the power grid and the algebraic equation reflecting the connection relationship of the power grid and the steady state process. The detailed blocks are as follows.

2.1) The physical side differential equation of the cyber-physical power systems under distributed control is established as follows.

{dot over (x)} ^(p)(t)=f ^(p)(x ^(p)(t),y ^(p)(t),u ^(p)(t))  (41)

where f^(p) denotes the dynamic equation of the power systems, x^(p)(t) denotes the vector of the physical state variables of the cyber-physical power systems at time point t and satisfies x^(p)(t)=[x₁ ^(pT)(t), x₂ ^(pT)(t), . . . , x_(N) ^(pT)]^(T), the superscript T denotes transpose, {dot over (x)}^(p)(t) denotes the differential of x^(p)(t) with respect to time, y^(p)(t) denotes the vector of the physical algebraic variables at time point t (since the power flow equations associated with physical algebraic variables require to be calculated uniformly throughout the grid, y^(p)(t) is not listed in nodes), and u^(p)(t) denotes the vector of the physical control variables at time point t and satisfies u^(p)(t)=[u₁ ^(pT)(t), u₂ ^(pT)(t), . . . , u_(N) ^(pT)(t)]^(T). 2.2) The physical side algebraic equation of the cyber-physical power systems under distributed control is established as follows.

0=g ^(p)(x ^(p)(t),y ^(p)(t))  (42)

where g^(p) denotes the stability power flow equation of the power systems.

2.3) The physical side algebraic equation with time-delay of the cyber-physical power systems under distributed control is established as follows.

$\begin{matrix} \left\{ \begin{matrix} {{0 = {{{g^{p}\left( {{x^{p}\left( {t - \tau_{m,i}} \right)},{y^{p}\left( {t - \tau_{m,i}} \right)}} \right)}\mspace{14mu} i} = 1}},2,\ldots\mspace{11mu},N} \\ {{0 = {{{g^{p}\left( {{x^{p}\left( {t - \tau_{e,i}} \right)},{y^{p}\left( {t - \tau_{e,i}} \right)}} \right)}\mspace{14mu} i} = 1}},2,\ldots\mspace{11mu},N} \end{matrix} \right. & (43) \end{matrix}$

2.4) The simultaneous equations of formulas (5)-(7) to establishing the physical power grid model of the cyber-physical power systems under distributed control are established.

3) The unified differential algebraic equation model of the cyber-physical power systems under distributed control is established.

The simultaneous equations of the formula (4) of the cyber side model in the block 1) and the formulas (5)-(7) of the physical power grid model in the block 2) are established, to obtain the unified differential algebraic equation model of the cyber-physical power systems under distributed control as follows.

$\begin{matrix} \left\{ \begin{matrix} {{{\overset{\cdot}{x}}^{p}(t)} = {f^{p}\left( {{x^{p}(t)},{y^{p}(t)},{u^{p}(t)}} \right)}} \\ {0 = {g^{p}\left( {{x^{p}(t)},{y^{p}(t)}} \right)}} \\ {{0 = {{{g^{p}\left( {{x^{p}\left( {t - \tau_{m,i}} \right)},{y^{p}\left( {t - \tau_{m,i}} \right)}} \right)}\mspace{14mu} i} = 1}},2,\ldots\mspace{11mu},N} \\ {{0 = {{{g^{p}\left( {{x^{p}\left( {t - \tau_{e,i}} \right)},{y^{p}\left( {t - \tau_{e,i}} \right)}} \right)}\mspace{14mu} i} = 1}},2,\ldots\mspace{11mu},N} \\ {{{{\overset{\cdot}{x}}_{i}^{c}(t)} = {{{h_{i}\ \left( {{z_{i}^{c}(t)},{x_{J \in N_{i}}^{c}\ \left( {t - \tau_{c,{ij}}} \right)}} \right)}\mspace{14mu} i} = 1}},2,\ldots\mspace{11mu},N} \\ {{0 = {{{u_{i}^{p}(t)} - {{C_{i} \cdot {x_{i}^{c}\left( {t - \tau_{e,i}} \right)}}\mspace{14mu} i}} = 1}},2,\ldots\mspace{11mu},N} \\ {0 = {{z_{i}^{c}(t)} - {dia{{g\left( {M_{x,i},M_{y,i}} \right)} \cdot \left\lbrack {{x_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T},{y_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T}} \right\rbrack^{T}}}}} \\ {{i = 1},2,\ldots\mspace{11mu},N} \end{matrix} \right. & (44) \end{matrix}$

4) The time-delay characteristic equation expression of the cyber-physical power systems under distributed control is established. The detailed blocks are as follows.

4.1) The equilibrium point of the cyber-physical power systems under distributed control is calculated, and the formula (8) is solved by making the differentials {dot over (x)}^(p)(t) and {dot over (x)}_(i) ^(c)(t) in the formula (8) both being 0, to obtain the vector of the equilibrium points (that is, a point at which the system may maintain stability without disturbance) of the cyber-physical power systems at time point t, namely, x(t)=[x ^(pT)(t), x ₁ ^(cT)(t), . . . , x _(i) ^(cT)(t), . . . , x _(N) ^(cT)(t)]^(T), where x ^(pT)(t) denotes the vector of the equilibrium points of the physical state variables at time point t (the vector includes the equilibrium point of each physical state variable), x _(i) ^(cT)(t) denotes the equilibrium point of the information state variable of the node i at time point t.

4.2) The linearization expansion is performed on the formula (8) of the unified differential algebraic equation model of the cyber-physical power systems under distributed control at the equilibrium point, and the following linearized differential algebraic equations may be obtained.

$\begin{matrix} \left\{ \begin{matrix} {{\Delta{{\overset{\cdot}{x}}^{p}(t)}} = {{A_{0}\Delta\;{x^{p}(t)}} + {B_{0}\Delta{y^{p}(t)}} + {\sum\limits_{i = 1}^{N}{E_{i}\Delta\;{u_{i}^{p}(t)}}}}} \\ {0 = {{C_{0}\Delta{x^{p}(t)}} + {D_{0}\Delta{y^{p}(t)}}}} \\ {{0 = {{{C_{e,i}\Delta\;{x^{p}\left( {t - \tau_{e,i}} \right)}} + {D_{e,i}\Delta\;{y^{p}\left( {t - \tau_{e,i}} \right)}\mspace{14mu} i}} = 1}},2,\ldots\mspace{11mu},N} \\ {{0 = {{{C_{m,i}\Delta\;{x^{p}\left( {t - \tau_{m,i}} \right)}} + {D_{m,i}\Delta\;{y^{p}\left( {t - \tau_{m,i}} \right)}\mspace{14mu} i}} = 1}},2,\ldots\mspace{11mu},N} \\ {{{\Delta{\overset{\cdot}{x}}_{i}^{c}} = {{{\sum\limits_{J \in N_{i}}{A_{ij}^{c}\Delta{x_{j}^{c}\left( {t - \tau_{c,{ij}}} \right)}}} + {B_{i}^{c}\Delta\;{z_{i}^{c}(t)}\mspace{14mu} i}} = 1}},2,\ldots\mspace{11mu},N} \\ {{{\Delta\; u_{i}^{p}} = {{C_{i}\Delta\;{x_{i}^{c}\left( {t - \tau_{e,i}} \right)}\mspace{14mu} i} = 1}},2,\ldots\mspace{11mu},N} \\ {{{\Delta z}_{i}^{c} = {{{K_{x,i} \cdot {x_{i}^{p}\left( {t - \tau_{m,i}} \right)}} + {{K_{y,i} \cdot {y_{i}^{p}\left( {t - \tau_{m,i}} \right)}}\mspace{14mu} i}} = 1}},2,\ldots\mspace{11mu},N} \end{matrix} \right. & (45) \end{matrix}$

where, Δ denotes a deviation value of the variable after performing the linearization expansion at the equilibrium point, A₀ denotes a Jacobi matrix obtained by taking derivative of f^(p) with respect to x^(p)(t), namely satisfies A₀=∂f^(p)/∂x^(p)(t), B₀ denotes a Jacobi matrix obtained by taking the derivative of f^(p) with respect to y^(p)(t), namely satisfies B₀=∂f^(p)/∂y^(p)(t), E_(i) denotes a Jacobi matrix obtained by taking the derivative of f^(p) with respect to u_(i) ^(p)(t), namely satisfies E_(i)=∂f^(p)/∂u_(i) ^(p)(t), i=1, . . . , N, C₀ denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to x^(p)(t), namely satisfies C₀=∂g^(p)/∂x^(p)(t), D₀ denotes the Jacobi matrix obtained by taking the derivative of g^(p) with respect to y_(p)(t), namely satisfies D₀=∂g^(p)/∂y^(p)(t), C_(e,i) denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to x^(p)(t−τ_(e,i)), namely satisfies C_(e,i)=∂g^(p)/∂x^(p)(t−τ_(e,i)), D_(e,i) denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to y^(p)(t−τ_(e,i)), namely satisfies D_(e,i)=∂g^(p)/∂y^(p)(t−τ_(e,i)), C_(m,i) denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to x^(p) (t−τ_(m,i)), namely satisfies C_(m,i)=∂g^(p)/∂x^(p)(t−r_(m,i)), D_(m,i), denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to y^(p)(t−τ_(m,i)), namely satisfies D_(m,i)=∂g^(p)/∂y^(p)(t−τ_(m,i)), A_(ij) ^(c) denotes a Jacobi matrix obtained by taking the derivative of h_(i) with respect to x_(j) ^(c)(t−τ_(c,ij)), namely satisfies A_(ij) ^(c)=∂h_(i)/∂x_(j) ^(c)(t−τ_(c,ij)), B_(i) ^(c) denotes a Jacobi matrix obtained by taking the derivative of h_(i) with respect to z_(i) ^(c)(t), namely satisfies B_(i) ^(c)=∂h_(i)/∂z_(i) ^(c)(t), C_(i) remains unchanged, K_(x,i) denotes a Jacobi matrix obtained by taking the derivative of diag(M_(x,i),M_(y,i)) with respect to x_(i) ^(p)(t−τ_(m,i)), namely satisfies K_(x,i)=∂diag(M_(x,i),M_(y,i))/∂x_(i) ^(p)(t−τ_(m,i)), and K_(y,i) denotes a Jacobi matrix obtained by taking the derivative of diag(M_(x,i),M_(y,i)) with respect to y_(i) ^(p)(t−τ_(m,i)), namely satisfies K_(y,i)=∂diag(M_(x,i),M_(y,i))/∂y_(i) ^(p)(t−τ_(m,i)).

4.3) Algebraic equations in the formula (9) of the linearized differential algebraic equations are eliminated, and the formula (9) is equivalently transformed into time-delay differential equations of the cyber-physical power systems under distributed control as follows.

$\begin{matrix} \left\{ \begin{matrix} {{\Delta{\overset{\cdot}{x}(t)}} = {{\begin{bmatrix} {A_{0} - {B_{0}{D_{0}}^{- 1}C_{0}}} & 0 \\ 0 & 0 \end{bmatrix}\Delta\;{x(t)}} + {\sum\limits_{i = 1}^{N}{{\begin{bmatrix} 0 & {E_{i}C_{i}} \\ 0 & 0 \end{bmatrix} \cdot \Delta}\;{x\left( {t - \tau_{e,i}} \right)}}} +}} \\ {{\sum\limits_{i = 1}^{N}{\sum\limits_{j \in N_{i}}\ {{\begin{bmatrix} 0 & 0 \\ 0 & A_{ij}^{c} \end{bmatrix} \cdot \Delta}\;{x\left( {t - \tau_{c,{ij}}} \right)}}}} +} \\ {\sum\limits_{i = 1}^{N}\ {{\begin{bmatrix} 0 & 0 \\ {{B_{i}^{c}K_{x,i}} - {B_{i}^{c}K_{y,i}{D_{m,i}}^{- 1}C_{m,i}}} & 0 \end{bmatrix} \cdot \Delta}\;{x\left( {t - \tau_{m,i}} \right)}}} \end{matrix} \right. & (46) \end{matrix}$

where the superscript ⁻¹ of the matrix denotes an inverse matrix of this matrix, x denotes the composite vector of the physical state variables and the information state variables, and satisfies x=[x^(pT), x₁ ^(cT), . . . , x_(i) ^(cT), . . . x_(N) ^(cT)]^(T).

4.4) The following time-delay characteristic equation of the cyber-physical power systems under distributed control is established.

$\begin{matrix} {{\det\begin{pmatrix} {{s \cdot I} - {\overset{\sim}{A}}_{0} - {\sum\limits_{i = 1}^{N_{c}}\;{T_{m,i} \cdot e^{{- s}\;\tau_{m,i}}}}} \\ {{- {\sum\limits_{i = 1}^{N_{c}}\;{\sum\limits_{j \in N_{i}}\;{T_{ij} \cdot e^{{- s}\;\tau_{c,{ij}}}}}}} - {\sum\limits_{i = 1}^{N_{c}}\;{T_{e,i} \cdot e^{{- s}\;\tau_{e,i}}}}} \end{pmatrix}} = 0} & (47) \end{matrix}$

where det denotes the determinant, s denotes the characteristic root, e denotes the natural exponential, and Ã₀, T_(ij), T_(m,i), T_(e,i) respectively satisfy:

$\begin{matrix} \left\{ \begin{matrix} {{{\overset{\sim}{A}}_{0} = \begin{bmatrix} {A_{0} - {B_{0}{D_{0}}^{- 1}C_{0}}} & 0 \\ 0 & 0 \end{bmatrix}},{T_{ij} = \begin{bmatrix} 0 & 0 \\ 0 & A_{ij}^{c} \end{bmatrix}}} \\ {{T_{m,i} = \begin{bmatrix} 0 & 0 \\ {{B_{i}^{c}K_{x,i}} - {B_{i}^{c}K_{y,i}{D_{m,i}}^{- 1}C_{m,i}}} & 0 \end{bmatrix}},{T_{e,i} = \begin{bmatrix} 0 & {E_{i}C_{i}} \\ 0 & 0 \end{bmatrix}}} \end{matrix} \right. & (48) \end{matrix}$

5) The time-delay of each node is determined and the stability criterion for the time-delay of the cyber-physical power systems under distributed control is performed. The detailed blocks are as follows.

5.1) The time-delay for measuring of each node τ_(m,i), i=1, 2 . . . , N, the time-delay for controlling of each node τ_(e,i), i=1, 2 . . . , N and the time-delay τ_(c,ij) for communication between each node i=1, 2, . . . , N and neighbor nodes j∈N_(i) of the node are set, then the formula (11) of the time-delay characteristic equation is solved, to obtain all characteristic roots of the equation. The characteristic root having the maximum real part is denoted as λ₁, and the real part of λ₁ is Re(λ₁).

5.2) It is determined whether Re(λ₁) is greater than 0. In response to Re(λ₁) being greater than 0, the cyber-physical power systems under distributed control is determined to be unstable. In response to Re(λ₁) being less than or equal to 0, the cyber-physical power systems under distributed control is determined to be stable. The stability criterion for the time-delay completes. 

What is claimed is:
 1. A stability criterion for time-delay of cyber-physical power systems under distributed control, comprising: establishing an cyber side model of the cyber-physical power systems under distributed control and a physical power grid model of the cyber-physical power systems under distributed control respectively; then establishing simultaneous equations of the cyber side model and the physical power grid model to establish an unified differential algebraic equation model of the cyber-physical power systems under distributed control, so as to obtain a time-delay characteristic equation expression of the cyber-physical power systems under distributed control; finally determining a time-delay of each node, and solving the time-delay characteristic equation expression to obtain a maximum characteristic root, and performing the stability criterion for the time-delay of the cyber-physical power systems under distributed control according to a real part of the maximum characteristic root.
 2. The stability criterion for time-delay of cyber-physical power systems under distributed control according to claim 1, comprising: 1) establishing the cyber side model of the cyber-physical power systems under distributed control; 1.1) establishing an equation for a measuring link of each node of the cyber-physical power systems under distributed control; wherein, the equation for the measuring link of a node i of the cyber-physical power systems under distributed control is as follows: $\begin{matrix} {{z_{i}^{c}(t)} = {\begin{bmatrix} M_{x,i} & 0 \\ 0 & M_{y,i} \end{bmatrix} \cdot \begin{bmatrix} {x_{i}^{p}\left( {t - \tau_{m,i}} \right)} \\ {y_{i}^{p}\left( {t - \tau_{m,i}} \right)} \end{bmatrix}}} & (1) \end{matrix}$ where a superscript c denotes an cyber side of the cyber-physical power systems, a superscript p denotes a physical side of the cyber-physical power systems, x_(i) ^(p)(t−τ_(m,i)) denotes a physical state variable of the node i at time point t−τ_(m,i), y_(i) ^(p)(t−τ_(m,i)) denotes a physical algebraic variable of the node i at time point t−τ_(m,i), z_(i) ^(c)(t) denotes an information measuring variable of the node i at time point t, τ_(m,i) denotes a time-delay for measuring of the node i, M_(x,i) denotes an incidence matrix from the physical state variable of the node i to the information measuring variable of the node i, and M_(y,i) denotes an incidence matrix from the physical algebraic variable of the node i to the information measuring variable of the node i; 1.2) establishing an equation for a communication iterative calculating link of each node of the cyber-physical power systems under distributed control; wherein, the equation for the communication iterative calculating link of the node i of the cyber-physical power systems under distributed control is as follows: {dot over (x)} _(i) ^(c)(t)=h _(i)(z _(i) ^(c)(t),x _(j∈N) _(i) ^(c)(t−τ _(c,ij)))  (2) where x_(i) ^(c)(t) denotes an information state variable of the node i at time point t, {dot over (x)}_(i) ^(c)(t) denotes a differential of x_(i) ^(c)(t) with respect to time, N_(i) denotes a set of reference numbers of neighbor nodes communicating with the node i each other, x_(j∈N) _(i) ^(c)(t) denotes an information state variable of the neighbor node of the node i at time point t, τ_(c,ij) denotes a time-delay for communication between the node i and a neighbor node j of the node i, and h_(i) denotes an iterative correction function of the node i under distributed control; 1.3) establishing an equation for a control link of each node of the cyber-physical power systems under distributed control; wherein, the equation for the control link of the node i of the cyber-physical power systems under distributed control is as follows: u _(i) ^(p)(t)=C _(i) ·x _(i) ^(c)(t−τ _(e,i))  (3) where C_(i) denotes an incidence matrix from the information state variable of the node i to a physical control variable of the node i, and u_(i) ^(p)(t) denotes the physical control variable of the node i, and τ_(e,i) denotes a time-delay for controlling of the node i; 1.4) establishing simultaneous equations of the equations for the measuring link, the equations for the communication iteration calculating link and the equations for the control link of all nodes together, and establishing the cyber side model of the cyber-physical power systems under distributed control as follows; $\begin{matrix} \left\{ \begin{matrix} {{{\overset{\cdot}{x}}_{i}^{c} = {{{h_{i}\ \left( {{z_{i}^{c}(t)},{x_{j \in N_{i}}^{c}\left( {t - \tau_{c,{ij}}} \right)}} \right)}\mspace{9mu} i} = 1}},{2\mspace{14mu}\ldots}\mspace{11mu},N} \\ {{0 = {{u_{i}^{p} - {{C_{i} \cdot {x_{i}^{c}\left( {t - \tau_{e,i}} \right)}}\mspace{20mu} i}} = 1}},{2\mspace{14mu}\ldots}\mspace{11mu},N} \\ {0 = {{z_{i}^{c}(t)} - {dia{{g\left( {M_{x,i},M_{y,i}} \right)} \cdot \left\lbrack {{x_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T},{y_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T}} \right\rbrack^{T}}}}} \\ {{i = 1},{2\mspace{14mu}\ldots}\mspace{11mu},N} \end{matrix} \right. & (4) \end{matrix}$ where diag( ) denotes a partitioned diagonal matrix, and N denotes the number of nodes for the distributed control; 2) establishing the physical power grid model of the cyber-physical power systems under distributed control; 2.1) establishing a physical side differential equation of the cyber-physical power systems under distributed control as follows; {dot over (x)} ^(p)(t)=f ^(p)(x ^(p)(t),y ^(p)(t),u ^(p)(t))  (5) where f^(p) denotes a dynamic equation of the power systems, x^(p) (t) denotes a vector of the physical state variables of the cyber-physical power systems at time point t and satisfies x^(p)(t)=[x₁ ^(pT)(t), x₂ ^(pT)(t), . . . , x_(N) ^(pT)(t)]^(T), a superscript T denotes transpose, {dot over (x)}^(p)(t) denotes a differential of x^(p) (t) with respect to time, y^(p) (t) denotes a vector of the physical algebraic variables at time point t, and u^(p)(t) denotes a vector of the physical control variables at time point t and satisfies u^(p)(t)=[u₁ ^(pT)(t), u₂ ^(pT)(t), . . . , u_(N) ^(pT)(t)]^(T); 2.2) establishing a physical side algebraic equation of the cyber-physical power systems under distributed control as follows; 0=g ^(p)(x ^(p)(t),y ^(p)(t))  (6) where g^(p) denotes a stability power flow equation of the power systems; 2.3) establishing a physical side algebraic equation with time-delay of the cyber-physical power systems under distributed control as follows; $\begin{matrix} \left\{ \begin{matrix} {{0 = {{{g^{p}\left( {{x^{p}\left( {t - \tau_{m,i}} \right)},{y^{p}\left( {t - \tau_{m,i}} \right)}} \right)}\mspace{14mu} i} = 1}},2,\ldots\mspace{11mu},N} \\ {{0 = {{{g^{p}\left( {{x^{p}\left( {t - \tau_{e,i}} \right)},{y^{p}\left( {t - \tau_{e,i}} \right)}} \right)}\mspace{14mu} i} = 1}},2,\ldots\mspace{11mu},N} \end{matrix} \right. & (7) \end{matrix}$ 2.4) establishing simultaneous equations of formulas (5)-(7) to establishing the physical power grid model of the cyber-physical power systems under distributed control; 3) establishing simultaneous equations of the formula (4) of the cyber side model in the block 1) and the formulas (5)-(7) of the physical power grid model in the block 2), to obtain the unified differential algebraic equation model of the cyber-physical power systems under distributed control as follows; $\begin{matrix} \left\{ \begin{matrix} {{{\overset{\cdot}{x}}^{p}(t)} = {f^{p}\left( {{x^{p}(t)},{y^{p}(t)},{u^{p}(t)}} \right)}} \\ {0 = {g^{p}\left( {{x^{p}(t)},{y^{p}(t)}} \right)}} \\ {{0 = {{{g^{p}\left( {{x^{p}\left( {t - \tau_{m,i}} \right)},{y^{p}\left( {t - \tau_{m,i}} \right)}} \right)}\mspace{14mu} i} = 1}},2,\ldots\mspace{11mu},N} \\ {{0 = {{{g^{p}\left( {{x^{p}\left( {t - \tau_{e,i}} \right)},{y^{p}\left( {t - \tau_{e,i}} \right)}} \right)}\mspace{9mu} i} = 1}},2,\ldots\mspace{11mu},N} \\ {{{{\overset{\cdot}{x}}_{i}^{c}(t)} = {{{h_{i}\ \left( {{z_{i}^{c}(t)},{x_{j \in N_{i}}^{c}\left( {t - \tau_{c,{ij}}} \right)}} \right)}\mspace{14mu} i} = 1}},{2\mspace{14mu}\ldots}\mspace{11mu},N} \\ {{0 = {{{u_{i}^{p}(t)} - {{C_{i} \cdot {x_{i}^{c}\left( {t - \tau_{e,i}} \right)}}\mspace{14mu} i}} = 1}},{2\mspace{14mu}\ldots}\mspace{11mu},N} \\ {0 = {{z_{i}^{c}(t)} - {dia{{g\left( {M_{x,i},M_{y,i}} \right)} \cdot \left\lbrack {{x_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T},{y_{i}^{p}\left( {t - \tau_{m,i}} \right)}^{T}} \right\rbrack^{T}}}}} \\ {{i = 1},{2\mspace{14mu}\ldots}\mspace{11mu},N} \end{matrix} \right. & (8) \end{matrix}$ 4) establishing the time-delay characteristic equation expression of the cyber-physical power systems under distributed control; 4.1) calculating an equilibrium point of the cyber-physical power systems under distributed control, and solving the formula (8) by making the differentials {dot over (x)}^(p)(t) and {dot over (x)}_(i) ^(c)(t) in the formula (8) both being 0, to obtain a vector of the equilibrium points of the cyber-physical power systems at time point t, namely, x(t)=[x ^(pT)(t), x ₁ ^(cT)(t), . . . , x _(i) ^(cT)(t), . . . , x _(N) ^(cT)(t)]^(T), x ^(pT)(t) denotes a vector of the equilibrium points of the physical state variables at time point t, x _(i) ^(cT)(t) denotes the equilibrium point of the information state variable of the node i at time point t; 4.2) performing linearization expansion on the formula (8) of the unified differential algebraic equation model of the cyber-physical power systems under distributed control at the equilibrium point, and obtaining the following linearized differential algebraic equations; $\begin{matrix} \left\{ \begin{matrix} {{\Delta{{\overset{\cdot}{x}}^{p}(t)}} = {{A_{0}\Delta\;{x^{p}(t)}} + {B_{0}\Delta{y^{p}(t)}} + {\sum\limits_{i = 1}^{N}{E_{i}\Delta\;{u_{i}^{p}(t)}}}}} \\ {0 = {{C_{0}\Delta{x^{p}(t)}} + {D_{0}\Delta{y^{p}(t)}}}} \\ {{0 = {{{C_{e,i}\Delta{x^{p}\left( {t - \tau_{e,i}} \right)}} + {D_{e,i}\Delta\;{y^{p}\left( {t - \tau_{e,i}} \right)}\mspace{20mu} i}} = 1}},\ldots\mspace{11mu},N} \\ {{0 = {{{C_{m,i}\Delta\;{x^{p}\left( {t - \tau_{m,i}} \right)}} + {D_{m,i}\Delta\;{y^{p}\left( {t - \tau_{m,i}} \right)}\mspace{14mu} i}} = 1}},\ldots\mspace{11mu},\ N} \\ {{{\Delta{\overset{\cdot}{x}}_{i}^{c}} = {{{\sum\limits_{j \in N_{i}}{A_{ij}^{c}\Delta{x_{j}^{c}\left( {t - \tau_{c,{ij}}} \right)}}} + {B_{i}^{c}\Delta\;{z_{i}^{c}(t)}\mspace{14mu} i}} = 1}},2,\ldots\mspace{11mu},N} \\ {{{\Delta\; u_{i}^{p}} = {{C_{i}\Delta\;{x_{i}^{c}\left( {t - \tau_{e,i}} \right)}\mspace{14mu} i} = 1}},2,\ldots\mspace{11mu},N} \\ {{{\Delta z}_{i}^{c} = {{{K_{x,i} \cdot {x_{i}^{p}\left( {t - \tau_{m,i}} \right)}} + {{K_{y,i} \cdot {y_{i}^{p}\left( {t - \tau_{m,i}} \right)}}\mspace{14mu} i}} = 1}},2,\ldots\mspace{11mu},N} \end{matrix} \right. & (9) \end{matrix}$ where, Δ denotes a deviation value of the variable after performing the linearization expansion at the equilibrium point, A₀ denotes a Jacobi matrix obtained by taking derivative of f^(p) with respect to x^(p)(t), namely satisfies A₀=∂f^(p)/∂x^(p)(t), B₀ denotes a Jacobi matrix obtained by taking the derivative of f^(p) with respect to y^(p)(t), namely satisfies B₀=∂f^(p)/∂y^(p)(t), E_(i) denotes a Jacobi matrix obtained by taking the derivative of f^(p) with respect to u_(i) ^(p)(t), namely satisfies E_(i)=∂f^(p)/∂u_(i) ^(p)(t), i=1, . . . , N, C₀ denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to x^(p)(t), namely satisfies C₀=∂g^(p)/∂x^(p)(t), D₀ denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to y_(p)(t), namely satisfies D₀=∂g^(p)/∂y^(p)(t), C_(e,i) denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to x^(p)(t−τ_(e,i)), namely satisfies C_(e,i)=∂g^(p)/∂x^(p)(t−τ_(e,i)), D_(e,i) denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to y^(p)(t−τ_(e,i)), namely satisfies D_(e,i)=∂g^(p)/∂y^(p)(t−τ_(e,i)), C_(m,i) denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to x^(p) (t−τ_(m,i)), namely satisfies C_(m,i)=∂g^(p)/∂x^(p)(t−τ_(m,i)), D_(m,i), denotes a Jacobi matrix obtained by taking the derivative of g^(p) with respect to y^(p)(t−τ_(m,i)), namely satisfies D_(m,i)=∂g^(p)/∂y^(p)(t−τ_(m,i)), A_(ij) ^(c) denotes a Jacobi matrix obtained by taking the derivative of h_(i) with respect to x_(j) ^(c)(t−τ_(c,ij)), namely satisfies A_(ij) ^(c)=∂h_(i)/∂x_(j) ^(c)(t−τ_(c,ij)), B_(i) ^(c) denotes a Jacobi matrix obtained by taking the derivative of h_(i) with respect to z_(i) ^(c)(t), namely satisfies B_(i) ^(c)=∂h_(i)/∂z_(i) ^(c)(t), C_(i) remains unchanged, K_(x,i) denotes a Jacobi matrix obtained by taking the derivative of diag(M_(x,i),M_(y,i)) with respect to x_(i) ^(p)(t−τ_(m,i)), namely satisfies K_(x,i)=∂diag(M_(x,i),M_(y,i))/∂x_(i) ^(p)(t−τ_(m,i)), and K_(y,i) denotes a Jacobi matrix obtained by taking the derivative of diag(M_(x,i),M_(y,i)) with respect to y_(i) ^(p)(t−τ_(m,i)), namely satisfies K_(y,i)=∂diag(M_(x,i),M_(y,i))/∂y_(i) ^(p)(t−τ_(m,i)); 4.3) eliminating algebraic equations in the formula (9) of the linearized differential algebraic equations, and equivalently transforming the formula (9) into time-delay differential equations of the cyber-physical power systems under distributed control as follows; $\begin{matrix} \left\{ \begin{matrix} {{\Delta{\overset{\cdot}{x}(t)}} = {{\begin{bmatrix} {A_{0} - {B_{0}{D_{0}}^{- 1}C_{0}}} & 0 \\ 0 & 0 \end{bmatrix}\Delta\;{x(t)}} + {\sum\limits_{i = 1}^{N}{{\begin{bmatrix} 0 & {E_{i}C_{i}} \\ 0 & 0 \end{bmatrix} \cdot \Delta}\;{x\left( {t - \tau_{e,i}} \right)}}} +}} \\ {{\underset{i = 1}{\overset{N}{\;\sum}}{\underset{j \in N_{i}}{\;\sum}{{\begin{bmatrix} 0 & 0 \\ 0 & A_{ij}^{c} \end{bmatrix} \cdot \Delta}\;{x\left( {t - \tau_{c,{ij}}} \right)}}}} +} \\ {\underset{i = 1}{\overset{N}{\;\sum}}{{\begin{bmatrix} 0 & 0 \\ {{B_{i}^{c}K_{x,i}} - {B_{i}^{c}K_{y,i}{D_{m,i}}^{- 1}C_{m,i}}} & 0 \end{bmatrix} \cdot \Delta}\;{x\left( {t - \tau_{m,i}} \right)}}} \end{matrix} \right. & (10) \end{matrix}$ where a superscript ⁻¹ of a matrix denotes an inverse matrix of this matrix, x denotes a composite vector of the physical state variables and the information state variables, and satisfies x=[x^(pT), x₁ ^(cT), . . . , x_(i) ^(cT), . . . x_(N) ^(cT)]^(T); 4.4) establishing the following time-delay characteristic equation of the cyber-physical power systems under distributed control; $\begin{matrix} {{\det\begin{bmatrix} {{s \cdot I} - {\overset{\sim}{A}}_{0} - {\sum\limits_{i = 1}^{N_{c}}\;{T_{m,i} \cdot e^{{- s}\;\tau_{m,i}}}}} \\ {{- {\sum\limits_{i = 1}^{N_{c}}\;{\sum\limits_{j \in N_{i}}{T_{ij} \cdot e^{{- s}\;\tau_{c,{ij}}}}}}} - {\sum\limits_{i = 1}^{N_{c}}\;{T_{e,i} \cdot e^{{- s}\;\tau_{e,i}}}}} \end{bmatrix}} = 0} & (11) \end{matrix}$ where det denotes a determinant, s denotes a characteristic root, e denotes a natural exponential, and Ã₀, T_(ij), T_(m,i), T_(e,i) respectively satisfy: $\begin{matrix} \left\{ \begin{matrix} {{{\overset{\sim}{A}}_{0} = \begin{bmatrix} {A_{0} - {B_{0}{D_{0}}^{- 1}C_{0}}} & 0 \\ 0 & 0 \end{bmatrix}},{T_{ij} = \begin{bmatrix} 0 & 0 \\ 0 & A_{ij}^{c} \end{bmatrix}}} \\ {{T_{m,i} = \begin{bmatrix} 0 & 0 \\ {{B_{i}^{c}K_{x,i}} - {B_{i}^{c}K_{y,i}{D_{m,i}}^{- 1}C_{m,i}}} & 0 \end{bmatrix}},{T_{e,i} = \begin{bmatrix} 0 & {E_{i}C_{i}} \\ 0 & 0 \end{bmatrix}}} \end{matrix} \right. & (12) \end{matrix}$ 5) determining the time-delay of each node and performing the stability criterion for the time-delay of the cyber-physical power systems under distributed control; 5.1) setting the time-delay for measuring of each node τ_(m,i), i=1, 2 . . . , N, the time-delay for controlling of each node τ_(e,i), i=1, 2 . . . , N and the time-delay τ_(c,ij) for communication between each node i=1, 2, . . . , N and neighbor nodes j∈N_(i) of the node, then solving the formula (11) of the time-delay characteristic equation, to obtain all characteristic roots of the equation, denoting a characteristic root having the maximum real part as λ₁, the real part of λ₁ being Re(λ₁); 5.2) determining whether Re(λ₁) is greater than 0; in response to Re(λ₁) being greater than 0, determining the cyber-physical power systems under distributed control to be unstable; in response to Re(λ₁) being less than or equal to 0, determining the cyber-physical power systems under distributed control to be stable; completing the stability criterion for the time-delay. 